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The light scalar mesons below 1 GeV as tetraquark states are studied in the framework 
of the flux-tube model, the multi-body confinement instead of the additive two-body 
confinement is used. From the calculated results, we find that the light scalar mesons, cr, 
K could be well accommodated in the diquark-antidiquark tetraquark picture in the flux- 
tube model and they could be color confinement resonances. The mass of the first radial 
excited state of is f019 MeV, which is close to the mass of /o(980). Whereas 

ao(980) can not be fitted in this interpretation. 
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1. Introduction 

The charged k, a scalar meson, was recently observed by BES collaboration.tH The 
Breit-Wigner mass and the decay width are obtained to be 826 ± 49ttl MeV and 
449±156^8f MeV, and the pole position is determined to be (764±63l^4)-i(306± 
149l8f ) MeV/c2. They are in good agreement with those of the neutral k: the 
mass and the decay width are 878 ± 23l 55 MeV and 499 ± 53tlr MeV, respectively, 
observed by the BES and other collaborations.^ 

The understanding of scalar mesons, which have the same quantum numbers as 
the vacuum, is a crucial problem in low-energy quantum chromodynamics (QCD) 
since they could shed light on the chiral symmetry breaking mechanism and pre- 
sumably also on confinement in QCD. Although many properties of scalar mesons 
have been studied for decades, it is still a puzzle for the understanding of the in- 
ternal stru cture of the scalar mesons. Their masses do not fit into the quark model 
predictions.12141 The flavor structures of these light scalar mesons below 1 GeV, 



1 



February 21, 2012 1:15 WSPC/INSTRUCTION FILE kappaN 



2 C. R. Deng, J. L. Ping and F. Wang 

ao(980), /o(980), a and k, are still an open questions. In the qq configuration, the 
p-wave relative motion between q and q has to be invoked to account for the spin 
and parity of the scalar mesons. This leads to much higher masses for them. Another 
possible configuration for scalar mesons is tetrquark state. In tetraquark configu- 
ration, the light scalar mesons could be classified into an SU(3) flavor nonet if the 
diquark picture is used.^ 6 7 Their quark contents can be expressed as 

- [su][su] + [sd\[sd] ^ 
a = [ud[ud\.,}^ = , 

= [ud\ [ds] , K+ — [ds] [ud] , — [ud] [us] , — [us] [ud] ; 

+ r ir-Ji [su][sd] - [sd][su] _ , 
% = MM],a^ = -j= ,ao = [sd][su]. 

Jaffe et al. interpreted light scalar mesons as tetr aquark states w ith all the rela- 
tive orbital angular momenta assumed to be zeros! * * III * I Weinstein et al. 
described light sc a.lar niesqns as hadroni c molecular states due to strong meson- 
meson interactionl ^^ l ^'^ l ^^ l ^^ l ^^ l ^^ l ^^ l ^^ l The properties of some of these light 
scalar mesons were also studied in the qq picture. The spectrum of light 

scalar mesons below 1.0 GeV were studied in the qq picture by including instan- 
ton interaction.'^ Bhavyashri et al. studied the instanton-induced interaction in 

light meson spectrum on the basis of the phenomenological harmonic models for 

^94 

quarks. Vijande et al. studied the scalar meso ns in te rms of the mixing of a chiral 
nonet of tetraquarks with conventional qq states .^2326] 

A multi-quark state as a multi-body system is quite different from ordinary 
hadrons {qq meson and qqq baryon) because the multi-quark state has more color 
structures than that of ordinary hadrons. The color structures of a multi-quark state 
is no longer trivial and the properties of the multi-quark states may be sensitive to 
the hidden color structure. A tetraquark state, if its existence is confirmed, may pro- 
vide important information about the low-energy QCD interaction which is absent 
from the ordinary hadrons. Some autho rs had studied the tetraquark system with 
the three-body qqq and qqq interaction.^^iEHl The exotic hadrons were also stud- 
ied as multiquark states in the flux tube model in our previous work.129 30 These 
studies suggest that the multi-body conflnement, instead of the additive two-body 
confinement, might be more suitable in the quark model study of multiquark states. 
The newly updated experimental data might shed more light on the possibility of 
existence of tetraquark states and QCD interaction for multi-quark states. 

The aim of this paper is to study the properties of scalar mesons below 1 GeV and 
their quark contents in the flux tube model with multi-body conflnement potential. 
The powerful method for few-body systems with high precision, Gaussian expansion 
method (GEM),^-*^ is used here. The paper is organized as follows: in Sec. 2, the 
flux-tube model with multi-body interaction are introduced. A brief introduction 
of GEM and the construction of the wave functions of tetra-quark states are given 
in Sec. 3. The numerical results and discussions are presented in Sec. 4 and a brief 
summary is given in the last section. 
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2. Quark model and multi-body confinement potential 

Long term studies in the past several decades on hadrons indicate that ordinary 
hadrons {qq meson and qqq baryon) can be well described by QCD-inspired quark 
models. Low energy QCD phenomena are dominated by two well known quark 
correlations: confinement and chiral symmetry breaking. The perturbative, effective 
one gluon exchange, properties of QCD should also be included. Hence, the main 
ingredients of the quark model are: constituent quarks with a few hundred MeV 
effective mass, phenomenological confinement potential, effective Goldstone bosons 
and one-gluon exchanged between these constituent quarks. 

For ordinary hadrons, the color structures of them are unique and trivial, naive 
models based on two-body color confinement interactions proportional to the color 
charges (Ai • \j) can describe the properties of ordinary hadrons well. However, 
the structures of a multiquark state are abundant j '^^ l '^^ l '^ ^ which include important 
QCD informations that is absent from ordinary hadrons. There is no any theoretical 
reason to extend directly the two-body confinement in naive quark model to multi- 
quark system. Furthermore, the direct application of the two body confinement to 
multi-quark system induces many serious problems, such as ant i- confinement 
and color Van der Waals force. Many theoretical works has been done to try to 
amend those serious drawbacks. The string flip model for mult i - quark system was 
proposed by M. Oka to avoid the pathological Van der Waals force. ^^^iSlxhree quark 
confinement is explored by introducing strings which connect quarks according to 
a certain configuration rule. ^^^^^^ 

Recent lattice QCD studies I 35 | 36 | 37 | gj^Q-^^ thai the confinement of multi-quark 
states is multi-body interaction and is proportional to the minimum of the total 
length of strings which connect the quarks to f orni a mu ltiquark state. Based on 
these studies, a naive flux-tube or string model l2 9 | 30 | 32 | multi-body conflne- 
ment has been proposed for multiquark systems. The harmonic interaction approxi- 
mation, i.e., the total length of strings is replaced by the sum of the square of string 
lengths, is assumed to simplify the numerical calculation. 

The diqu ark-antidiq uark picture of tetraquark states has been discussed by sev- 
eral authors .^^EHISQllO] ^Yiq present work, the scalar mesons below 1 GeV are 
studied as diquark-antidiquark systems in the flux-tube model. The color structure 
is shown in Fig.l, where a solid dot represents a quark while a hollow dot represents 
an anti-quark, r.^ is quark's position, represents a junction where three strings 
(flux tubes) meet. A thin line connecting a quark and a junction represents a fun- 
damental representation, i.e. color triplet. A thick line connecting two junctions is 
for a color sextet or other representations, na mely a com pound string. The differ- 
ent types of string may have different stiffness.SlE^E31 Color couplings satisfying 
overall color singlet of tetra-quark are [[^^Isi^^lsli and [['Z'zleiwlelii the subscripts 
represent the dimensions of color representations. 

In the flux-tube model with quadratic confinement, the confinement potential, 
which is believed to be fiavor independent, of the tetraquark state has the following 
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form}^ 

V"" = k [(n - y,f + {V2 - yi)' + (ra - ya)' + [vi - ys)' + Kdivi - 72)'] , (1) 

where k is the stiffness of the string with the fundamental representation 3 which is 
determined by meson spectrum, kKd is the compound string stiffness. The compound 
string stiffness parameter k^I^ depends on the color representation, d, of the string, 

«d = ^, (2) 

where Cd is the eigenvalue of the Casimir operator associated with the SU{3) color 
representation d of the string. C3 = |, Ce = ^ and Cg = 3. 





Fig. 1. Diquark-antidiquark state. 

For given quark positions (i = 1, 4), we can fix the positions of the junctions 
Yi {i — 1,2) by minimizing the energy of the system. After fixing y^, a set of 
canonical coordinates Ri {i = 1,...,4) is introduced to simplify expressions of the 
potential, are read as. 



R3 



Ri = 



R2 



r2 + rs + r4). 



r2 - rs - r4), R4 

Therefore, the minimum of the confinement interaction has the following form. 



1 + Kd 



Rg 



(3) 



Taking into account potential energy shift, the confinement potential l^^j„ used in 
the present calculation has the following form 



V 



c 

rain 



(R?-A). 



(R^ 



A) 



1 



Hd 



-(R3-A) 



(4) 



Carlson and Pandharipande also considered similar flux tube energy shift which 
they assumed to be proportional to the number of quarks Obviously, the 

confinement potential V'-^ is a multi-body interaction rather than a two-body in- 
teraction. It should be emphasized here that our approach is different from that in 
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Iwasaki's workp^ where the four-body problem is simplified to two-body one by 
treating diquark as a antiquark and antidiquark as a quark. 

The other parts of the Hamiltonian are rest masses, kinet ic en ergies, one-gluon- 
exchange potential and Goldstone-boson-exchange potentials ,^211 



i>3 



(5) 



9ch 



A2 

X 



47r 12mimj — 



A3 



X = TT,K,ri. 

(6) 



where Tcm is the center-of-mass kinetic energy, Fi,Xi are flavor, color SU3 Gell- 
mann matrices. Y{x) is the standard Yukawa function an d all o ther symbols have 
their usual meaning. The (5- function should be regularized ^^SHZl 



47r njrl{ij) 



(7) 



where is the reduced mass of qi and qj and rn(/i) = fQ/ji. The effective scale- 
dependent strong coupling constant is given by^^ 



"0 



In 



(8) 



3. Wave functions and gaussian expansion method 

The total wave function of a tetra-quark state can be written as a sum of the 
following direct products of color, isospin and spatial-spin terms, 

^IMiJtMt = E ^li^i.cX 
l.s,i.c.L 

(9) 

Here / and Jt are total isospin and total angular momentum respectively, ^'sim^j^ 
) and Xciwi {XC2W2) are spin, flavor and color wave func- 
tions of, diquark (anti-diquark), respectively. [ ]'s denote Clebsh-Gordan coefficients 
coupling. The coefficient ^^"^^^ c l determined by diagonalizing the Hamiltonian, 
subscripts s, i, c, L represent all possible intermediate quantum numbers, therefore 
our calculations are multi-channel coupling calculations. The Jacobi coordinates of 
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the tetraquark are defined as 

r = ri - r2, R = ra - r4 

TOiTi + m2r2 r7i3i"3 + miVi miri + m2r2 + m^T^ + m4YA 
A - ■ ■ , KcM - ; ; ; 

mi + 1712 "I3 + m4 TOi + TO2 + TO3 + m4 

where particles 1 and 2 are two quarks and particle 3 and 4 are two anti-quarks. 
L, li and I2 are the orbital angular momenta associated with coordinates X, r and 
R, respectively. The calculation is done in the center-of-mass coordinate system 
(RcM — 0). The tetra-quark state is an overall color singlet with well defined 
parity P — (— l^i+'a+i^ isospin / and total angular momentum Jt- For scalar 
mesons, we set the angular momentum L, li and I2 to be zero. 

For color part, the color singlet is constructed in the following two ways, x\ = 
3i2 ® 834, Xc = 612 ® 634, both "good" diquark and "bad" diquark are included. 
With respect to the flavor part, the flavor wave function read as 77/ = 7712 ® rj^^. 
Taking into account all degrees of freedom, the Pauli principle must be satisfied 
for each subsystem of identical quarks or anti-quarks. To obtain a reliable solution 
of few-body problem, a high precision method is indispensable. In this work, the 
GEM,SlI which has been proven to be rather powerful to solve few-body problem, 
is used to do the calculations. In GEM, three relative motion wave functions are 
written as, 

"Imax 

ni — 1 

i^hmA'^)= 51 c„,^„,i,i?'^e-'^"^^Vi,™,(R) (10) 



"2=1 



"3 = 1 

Gaussian size parameters are taken as the following geometric progression numbers 



4. Numerical results and discussions 

Now we turn to the tetraquark state calculation. The model parameters are fixed 
by reproducing the ordinary meson spectrum and are listed in Table 1. the meson 
spectrum can be reproduced very well. Because the flux-tube model is reduced to 
the ordinary quark model for qq system, the obtained meson spectra (from light 
to heavy) are similar to that of other work, e.g. Ref-TH Parts of the calculated 
meson spectrum are shown in Table 2. The experimental values are taken from 
PDG compilation.SHI 
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Table 1. The model parameters (set I). The masses of 7r,K,ri take the experimental values. 



m ma k fo Ao f^o A ao A,r Ak = A,, dp 

MeV MeV MeV fm-^ MeV fm fm^^ fm^^ fm^ - fm~i fm^^ 
313 520 213.3 30.85 0.187 0.113 0.5 4.25 4.2 5.2 15° 



Table 2. The meson spectrum (unit: MeV). 



Meson 


7r 


K 


P 


K* 







Cal. 


139 


502 


761 


897 


735 


1023 


Exp. 


139 


496 


770 


898 


780 


1020 


/(?2y(fm) 


0.57 


0.60 


1.05 


0.96 


1.02 


0.85 



The energies of scalar meson states can be obtained by solving the four-body 
Schrodinger equation 

= (12) 

with Rayleigh-Ritz variational principle. In GEM the calculated results are con- 
verged with nimax—(i, n2max = 6 and 713^02: = 6. Minimum and maximum ranges 
of the bases are 0.1 fm and 2.0 fm for coordinates r, R and X, respectively. 

Quark contents and the corresponding masses in the three different quark mod- 
els for the light scalar mesons as tetra-quark states are shown in Table 3, where 
n stands for a non-strange quark [u or d) while s stands for a strange quark, 
I and N denote total isospin and principal quantum number of total radial ex- 
citation, S*, L and J have their usual meaning. "Naive" stands for naive quark 
model, where only one-gluon-exchange potential is taken into account in addition 
to the additive two-body confinement.'^ "Chiral" stands for chiral quark model, 
where one-gluon-exchange and one-Goldstone-boson-exchange are included besides 
the additive two-body confinement. S3 The masses in the naive and chiral model are 
much higher (several hundreds MeV) than those in flux-tube model, the origin of 
this discrepancy mainly comes from the different type of confinement interaction, a 
two-body confinement potential is applied in the naive and chiral model, whereas a 
multi-body interaction confinement is used in the fiux-tube mode. Zou et al. studied 
scalar mesons in the quark model by introducing three-body confinement interac- 
tion. Their study also indicate that the multi-body confinement potential, inst ead 
of two-body interaction, should be applied in the study of multi-quark states.'^ 
The naive quark model gives the highest masses, due to the absence of Goldstone 
boson exchange, which induces additional attraction for the tetraquark system. 

In the framework of the flux tube model, it can be seen from Table 3 that the 
lowest masses of nnhn and nnns system are 587 MeV and 948 MeV, which are 
close to the masses of a and k. If the existence of cr and k is further confirmed, 
the tetraquark state is a possible i nterp retation. This interpretation is in agreement 
with many other studies.i ^^'*^ '^^■^''^^l'*''*^ ! Prelovsek et al. recently studied light scalar 
mesons a and k by lattice QCD simulation, they also found that a and k have sizable 
tetra-quark components nnnh and nnns, respectively.^^ In the case of /o(980) and 
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Table 3. Numerical results for three models (unit: MeV). 



Flavor 
IJP 


nnnn 
00+ 
01 So 


nnnn 
00+ 

l^So 


nnnn 
10+ 
O^So 


nnns 
O^So 


nsns 
00+ 
01 So 


nsns 
10+ 
01 So 


Naive 


938 


1431 


1431 


1216 


1456 


1456 


Chiral 


666 


1237 


1406 


1122 


1454 


1454 


Flux-tube 


587 


1019 


1210 


948 


1314 


1318 


Candidate 
mass 


a 

541 ± 39 51 


/o(980)? 
980 ±10^ 




K 

826±49t^9|l] 







Table 4. The parameter set II. 



k 




A 




MeV fm~2 


MeV fm 


fm2 




267 


30.0 


0.6 


4.09 



Table 5. Numerical results in flux-tube model (unit: MeV). 



Flavor 


nnnn 


nnnn 


nnnn 


nnns 


nsns 


nsns 


IJP 


00+ 


00+ 


10+ 




00+ 


10+ 




OiSo 


liSo 


OlSo 


OiSo 


OiSo 


OiSo 


Flux-tube 


531 


969 


1180 


908 


1270 


1275 



ao(980), many theoretical studies assumed them as tetra-quark states with quark 
content nsns with isospin 7 = and 7=1, respectively. The masses for nsns 
states are much higher than experimental values even in the flux tube model. Thus 
their main components seem to be not the nsns in the quark model. Instead, the 
mass of the radial excited state of a is 1019 MeV, which is close to the mass of 
/o(980). Taking /o(980) as nnnn state is consistent with Vijande's work.l^Sl The 
observed /o (980) — >■ KK process can be explained by nsns components mixing and 
the calculation is on going. Pelaez also suggested that three scalar mesons, a, k and 
/o(980), have dominant tetraquark component, whereas ao(980) might be a more 
complicated systems,'''^ which is also consistent with our results. 

In order to check the dependence of numerical results on model parameters, 
we make the same calculations on scalar mesons with another set of parameters 
which listed in Table 4 (the unchanged parameters are not listed) . The almost same 
meson spectrum is obtained. The results for tetraquark states are shown in Table 
5. Comparing Table 4 and 5, our results are rather stable against the variation of 
model parameters. 

5. Summary 

The scalar mesons below 1 GeV are studied as diquark-antidiquark system in the 
flux-tube model including multi-body confinement interaction on the basis of the 
lattice QCD calculations and other two models involving two-body confinement po- 
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tential. In the naive and chiral quark models, it is difRcult to explain the observed 
scalar mesons as qq or tetra-quark states. In the flux-tube model, rather low masses 
of tetraquark states are obtained. It is possible to assign the low- lying scalar mesons 
as tetraquark states, at least tetraquark is the main component. For example, a and 
K can be assigned as nnnn and nnns with = 0^ , if their existence is further con- 
firmed. To assign /o(980) and ao(980) as the tetraquark states qsqs is questionable 
since the theoretical masses are much higher than their experimental values. The 
interpretation of the ao(980) as a lowest qqqq state would give a much higher mass 
(about 200 MeV higher) than experiment data, which is still a mysterious system 
for the time being. However, in our calculation the mass of the radial excitation 
of the cr is close the mass of /o(980). The problem of this assignment, the small 
decay width of /o(980) — >■ KK, can be accounted for by the mixing of nshs com- 
ponent with nnnn. Our results are rather stable against the variation of the model 
parameters. 

At now the nature of scalar meson is still an open question. The interpretation 
of the scalar mesons as tetraquark states is a possibility. In fact, the scalar mesons 
should be the superpositions of qq, qqqq and other components in a Fock space 
expansion approach and the dominant one is determined by quark dynamics. The 
mixing between two-body and four-body configurations would require the knowledge 
of the quark-antiquark pair creation-annihilation interaction, which is calculating in 
our group by using a '^Pq model tentatively. The crucial test of the quark structure 
of the scalar mesons should be determined by the systematic study of these scalar 
mesons' decay pattern. Of course, the final answer could only be obtained through 
a combined effort of the further experimental measurement and theoretical efforts.. 

As mentioned before, the flavor structures of scalar mesons are complicated, 
even limiting to four-quark system. At the same time, their color structures are also 
complicated, such as3iX)3, 6iX)6, 1(8)1 and 8 ig) 8. The mixing of different color 
structures should be determined by the transition interaction between different color 
structures. Unfortunately, at present we do not have any reliable information about 
these transition interactions. 

The comparative studies of three models reported here indicate that the flux- 
tube model with a multi-body confinement potential on the basis of lattice QCD can 
give lower masses of scalar mesons as tetra-quark states than other two models with 
two-body confinement interaction. The multi-body interaction plays an important 
role in reducing the energy of the system. We therefore strongly suggest that a multi- 
body confinement potential, instead of two-body interaction, should be employed in 
the study of multi-quark system. In this way, multi-quark systems can form multi- 
body flux-tubes connecting quarks into a color singlet multi-quark states. The states 
can not decay into two colorful hadrons directly due to color confinement. They must 
transform into color singlet mesons before decaying. This decay mechanism is similar 
to compound nucleus formation and therefore should induce a resonance which is 
named as a "color confined, multi-quark resonance" state before.'^'-' It is different 
from all of those microscopic resonances discussed by S. Weinberg. Currently, an 
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urgent task is to obtain the information about the transition interaction between 
different Fock space configurations and different color structures which are crucial 
both in the structure and decay width calculations. 

We apologize to those authors who made contributions in this field and their 
works have not been mentioned here due to the space limitations. This work is 
supported partly by the National Science Foundation of China under Contract Nos: 
11047140, 11175088, 11035006 and the Ph.D Program Funds of Chongqing Jiaotong 
University. 
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